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Abstract 

We study the phase distribution around a vortex in uniform motion. We 
consider both the cases of neutral and charged superfluids. The motion of 
the vortex causes the density of the system to fluctuate. This in turn pro- 
duces a compensating current, thus ensuring current conservation locally. We 
explicitly calculate this current. 
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In recent years, there has been a great deal of interest (see Ref. 1 for a recent review) in 
the static and dynamic properties of vortices in superfluids and superconductors. A study 
of high temperature and other superconductors has led to the discovery of several new 
phenomena such as flux-lattice melting (2), quantum flux creep (3), anomalous sign change 
of the Hall co-efficient (4) and anomalous a.c. electromagnetic response(5). This has led to 
a renewed upsurge in efforts to understand the dynamics of moving vortices. Issues like the 
size and sign of the Magnus force (6) and the vortex inertial mass (7)have been intensely 
studied. Theoretical studies of these phenomena have been carried out, both within an 
effective Ginzburg-Landau theory (8) and microscopically (9) within the BCS theory. 

However, most studies of vortex dynamics have tended to assume the vortex motion 
to be adiabatic i.e. the instantaneous distribution of the superconducting order parameter 
around the vortex is assumed to be identical with that of a vortex at rest. In this paper we 
address the question of non-adiabatic distortions in the phase distribution around a moving 
vortex. Our calculations are performed using a Ginzburg-Landau phase-only functional. 
We solve the equations of motion to second order in the vortex velocity (u) and thus obtain 
corrections to the phase distribution. 

The vortex motion induces fluctuations in the instantaneous density distribution around 
the vortex. Current conservation therefore causes a compensating current to flow, thus 
conserving the density locally. We explicitly calculate the induced current, to second order 
in u for neutral as well as charged superfluids. The constant current contours are found to 
be distorted from their circular shape. However, there is no change in the energy, to second 
order in u, coming from the non-adiabatic distortions in the case of a neutral superfluid. 
For a charged superfluid the displacement current, associated with a time-dependent vector 
potential, contributes to a change in the energy. This contribution is suppressed by a factor 
of u 2 /c 2 and is therefore negligibly small. These effects, thus have no consequences either 
for the Magnus force or for the vortex mass. 

The superfluid is described by the action functional per unit length, S' = S/L = 
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Here 6 1 is the phase of the superconducting order parameter and A Q and A are the scalar 
and vector potentials associated with the electromagnetic field. This functional has been 
derived microscopically (10) and the co-efficients a.\ and a 2 can be related to appropriate 
polarizabilities of the underlying fermionic system. The size and sign of the co-efficient 7, 
which determines the Magnus force, remains controversial. For purposes of this paper, we 
treat the Ginzburg-Landau coefficients as phenomenological parameters. 

We will first consider the case of a neutral superfluid. Accordingly, we put e = in Eq. 
(1) and switch off the electromagnetic fields. We then have, 

L ns = ^1q2 _^2 { y &) 2 +1 q 

The Euler-Lagrange equation of motion is easily found to be 

ol\Q = a 2 V 2 6 (3). 

Notice that the term linear in 9 corresponds to a topological phase and makes no contribution 
to the equation of motion. 

For a uniformly moving vortex we have 0(f,t) = 6{f — ut). Substituting this in Eq. (3) 

we get 

a x u ■ V(« • V0) = a 2 V • (V0) (4) 
We seek a solution to Eq. (4) of the form 

6(f,t) = e^\r-ut) + 6 {l \r-ut)+e {2) {r-ut) + .... (5) 

where 6^ rn \f—ut) is of m th order in u. Here V9^\r') = ^ (where we define the co-ordinate 
r' —f—ut for convenience) is the phase distribution of a static vortex while the remaining 



terms are non-adiabatic phase distortions induced by the vortex motion. The quantiza- 
tion of the vorticity, which stems from the single-valuedness of the superconducting order 
parameter, ensures that for m ^ 9^ m \r') is non-singular and V^"^ is purely longitudinal. 

Substituting Eq. (5) in Eq. (4) and equating terms of the same order in u, we arrive at 
the result 

V' 2 #«(r') = (6a) 

and 

ai u ■ V(u ■ V'fl (0) (?)) = a 2 V 2 9 {2) (?) (66) 

Since V9^ is constrained to be purely longitudinal, Eq. (6a) implies that it is zero. We 
therefore set 9^ = 0. 

Fourier transforming Eq. (6b) with respect to r' we solve for 9^ to get 

9 { > = u ■ q 7a 

and find the corresponding current V9^ to be 

v7//i(2) 27r ai u-zxq 

v#V = q( u -o) — — — ( 7b ) 

ct 2 %q^ 

It is easy to see from Eq. (2) that this current, which is purely longitudinal, doesn't con- 
tribute to the action to second order in u. This implies that both the Magnus force and 
the vortex mass are unaffected by the existence of this current. To get a better idea of 
the induced current distortion, we Fourier transform V'6* q (the details are outlined in the 
Appendix) to get 

V9{?) = ^ [1 - cos 20'] (8) 
r 2a 2 

for the x-axis chosen to be directed along u. It is easily seen from Eq. (8) that the circular 
symmetry of a static vortex is destroyed by the extra current. 

We now turn our attention to the charged case. The Euler-Lagrange equations of motion 
are now given by 

ai §i ( ° ~ 7T Ao) = a ' V ' (W " Yc A) (9a) 



and 
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We choose the gauge V • A = 0. Then, on Fourier transforming Eq. (9b) we find 

*«> = TeT^fr (10) 

where the Thomas-Fermi screening length is given by 47rai(^) 2 = \ T 2 F . On Fourier trans- 
forming Eq. (9c) and making use of Eq. (10) we have 

7,* (-£-g)gM-W q (ftc/2e)W q 

W 2ec(g 2 A^ + l)(g 2 + AZ 2 -(M-g7c) 2 ) g 2 A| + 1 - (A L « • g/c) 2 1 ; 

where the penetration depth \l is given by 4na2(2e/hc) 2 = Xj 2 . In writing Eq. (11), we 

— * — * 

have made use of the relations A(f,t) = A(f — ut) and A (f,t) = A (f — ut). All Fourier 

— * 

transforms have been performed with respect to the co-ordinate r' = f — ut as before. 

We now substitute our result for A (Eq.10) in Eq. (9a) and make use of our gauge 
condition to arrive at 

ai u-qu-V% 2 f 2 ^ - = a 2 q- V'0 q (12). 
q 2 Y TF + 1 

We put 9 = # (0) + + # (2) + as before. We once again find to be zero and get 



e m = -^ a . f->*t £™ (13) . 

q a 2 q 4 q 2 \ T F + 1 

Notice that as e — >• and \ TF — > oo, Eq.(13) recovers the earlier result of Eq. (7a) for a 
neutral superfluid. 

Substituting Eq. (13) in Eq. (11) we find, to second order in u, the vector potential to 

be 



he 2-k z x q A| (u ■ q)' 

2e~i^q 2 \ 2 L + V + ^ 2 ~g 2 A| + l J 

As expected, the vector potential is purely transverse as the longitudinal terms exactly 
cancel, consistent with our gauge condition. 
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Making use of Eqs. (13) and (14) we find that the current density j(r') = [V9 — j^A] is 



given by j = j a + j b + j c where 



2n zx qq 2 X\ 

3a{q) = — 2 \2 i 1 15a 
iq l q z X L + 1 



is the cuurent density of a vortex at rest, 



, . 2irX 2 r zxqq 2 X 2 r ._ , 2 . 
corresponds to the "displacement current" of the moving flux and 

Jc( ^ = W^-^^Ti (15c) 

is the induced current necessary for local charge conservation. 

We now proceed to evaluate j(r'). Fourier transforming j a (q) we find (see Appendix for 
details) 

j a (P) = Y~K±(r' / Xl) (16) 

which is the usual result for a vortex at rest. Here K m (x) are the modified Bessel functions. 
In a similar manner, we evaluate jbij-') and j c (r f ), details of which have been provided in 
the Appendix. 
We find 

- - u 2 

J b (r>) = —[((2 + cos20')/i - cos207 2 )0' + sin20'(/ 1 + I 2 )f'] (17) 



where 



and 



Similarly, 



Here 



h = ±-[3K^/X L ) + y L K[(r'/X L )] (18a) 



r' 

2A| 



K 2 (r'/X L ) (186). 



2 \ 2 

= \^\ih - h) sm20'r' - (J 3 + h) cos20'0'] (19). 



h ~ (20a), 



and 

h - ^(1 - Jo(r'/0) " ^Mr'/O + J^W/O (206), 

and J m (%) are Bessel functions. The smallness of J 3 and J 4 stems from the small screen- 
ing length (\tf) associated with j c (see Eq.(15c)). This length is much smaller than the 
coherence length £, which is the coarse-graining scale beyond which a description in terms 
of a phase-only action functional is applicable. Notice also that both j b and j c are strongly 
suppressed by factors of u 2 /c 2 making their observation a difficult task. 

We now consider the drawbacks of our calculation. As pointed out by Aitchison et. 
al. (10), at T = 0, one has to include additional terms in the phase-only Lagrangian to 
ensure Galilean invariance. However, the solution of the corresponding non-linear problem 
is beyond the scope of our work. We have also ignored effects arising from the vortex core, 
whose proper inclusion would require a use of microscopic theory. 

Finally, we summarise the main results of this paper. We have calculated the non- 
adiabatic phase distortion of a moving vortex. Vortex motion induces time-dependent density 
fluctuations, which in turn gives rise to additional currents necessary to ensure local charge 
conservation. For the neutral superfluid, the extra supercurrent is purely longitudinal with 
observable consequences. In the charged case there are two distinct contributions. One 
coming from the "displacement current" (j b above) is purely transverse whereas the other 
(j c above) is longitudinal. These currents are however strongly suppressed by factors of 
u 2 /c 2 as well as efficient Coulomb screening in the case of j c , making their observation a 
difficult task. 
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APPENDIX: 

In this appendix, we outline the derivation of the results stated above for the current 
j(r'). We will first take up the case of a neutral superfluid. Then we have 

2nz x q 1hol\ u-zxq 

V# q = — ; 2 ~t~ q{u-q) — — (Al) 

iq 1 a 2 iq A 

where the first term on the R.H.S. is the usual current of a vortex at rest and the second 
term is the distortion induced by the vortex motion. Thus, the current in real space is given 
by 

V9{?) = J ^- exp (iq- f')V'0 q (A2). 
The first term in Eq. (Al) is easily Fourier transformed to yield the standard result 

V9(r% = ^ (A3). 

We find, on performing the angular integrals in Eq. (A2), the second term in the current to 
be 

2 

V '0(r') 6 = / dq[f'sm29(J 3 (qr') - J^qr')) + 0' cos 29 (J 3 (qr>) + J^qr')) (AA) 

4q; 2 J 

where 9 is the angle u makes with respect to r' . The integrals in Eq. (A4) can now be 
performed to yield 

V9(r% = cos 29^- (Ah). 

v 1 2a 2 r' v 1 

On adding the results of Eqs. (A3) and (A5) and choosing u to lie along the x-axis we finally 
arrive at Eq.(8) above. 

We now consider the charged case. In this case, the current j(q) is the sum of three 
different contributions:^, jb and j c . The corresponding Fourier transforms j a (q), jb(o) an d 
j c (q) have been obtained earlier (Eq. (15)). We will evaluate each piece seperately. 

is the usual current associated with a static vortex. On performing the angular integration 
we get 

Ur') =$1 dq q ^^ = f-Ktf/Xz.) (A7) 
J q + \ L Al 
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We now consider jb(r'). 

dq - a 2ttA| z x qq 2 X\ 



On performing the angular integral we arrive at Eq. (17) above where 



and 



I\ can be rewritten as 



g 4 Aj Uqr) 



r _ rlA U / Jl{qr,) - z^iAr^VMi m 10 ) 

h ~ 2 9A L i % 2 A| + 1 " 2 9A L [ A| J (A1U) 

which reduces to the result of Eq. (18a). I 2 is directly evaluated to yield the result stated 
earlier in Eq. (18b). 

Finally we consider j c {r'). 

f dq -,s2n\ 2 , u-zxq , A ~~\ 

Jc(r } = J (2^ eXP {iq • r } ' ^T^Ti {All) 

On performing the angular integration in Eq. (All) we arrive at the result of Eq. (19) 
where 

g 2 Ji(gr') 



and 



7 3 is easily evaluated to yield 



q 2 J 3 (qr') 
q 2 X 2 TF + l 



dq^+ {Am). 



h = ^-K^t'/Xtf) (A13). 

The Thomas- Fermi screening length \ T f "C £ where £ is the superconducting coherence 
length beyond which our coarse-grained picture based on a phase-only functional is valid. 
Thus Ki(r' j Xtf) ~ ex P (—t'/Xtf) — > in this (r > £)regime. 
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To evaluate I4 we make use of the relation J m _i(x) — J m+ i(x) = 2J' m (x) which is obeyed 
by Bessel functions. Using this relation, we find that I4 can be rewritten as 

f q 2 J 1 (qr') d [ q[J (qr') - 2J[(qr')} 

h = l dq q ^ TF + 1 ~ 2 d?j dq — — {AU) - 

This expression can be further simplified to yield 

We ignore the first two terms on the right hand side of Eq. (A15) as they are exponentially 
small and in view of the small screening length, approximate the last integral by 

hmrk K f dqMqr ' } = ^ [1 - wm (m - 

Substituting these approximations in Eq. (A15) we finally arrive at the result of Eq. (20b). 
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